Optimizing the search for resources by sharing information: Mongolian gazelles as a 

case study. 
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We investigate the relationship between communication and search efficiency in a biological con- 
text by proposing a model of Brownian searchers with long-range pairwise interaction. After a 
general study of the properties of the model, we show an application to the particular case of acous- 
tic communication among Mongolian gazelle, for which data are available, searching for good habitat 
areas. Using Monte Carlo simulations and density equations, our results point out that the search 
is optimal (i.e. the mean first hitting time among searchers is minimum) at intermediate scales of 
communication, showing that both an excess and a lack of information may worsen it. 

PACS numbers: 05.40.-a, 05.40.Fb, 87.23.Cc 



Many living organisms, including bacteria insects, 
and mammals 2, 3] communicate for a variety of reasons 
including to facilitate social cohesion [3, [![, to defend 
against predators @, to maintain territories SB, and 
to share information on resource locations [HT HI- Such 
strategies often lead to spatial group formation [10j and 
may be mediated by a range of mechanisms, including 
vocal communication 11] . Group formation strategies, 



such as individuals who have located resource rich areas 
attempting to attract conspecifics to help reduce their 
own vulnerability to predators, may also have important 
incidental benefits such as improving average foraging ef- 
ficiency. This kind of strategy might be favored because it 
is a "win-win" situation; the caller benefits from reduced 
predation risk, while individuals responding to the call 
learn the locations of high resource areas. 

A variety of mammalian species are known to com- 
municate acoustically over distances up to several kilo- 
meters 0], but while group formation via vocalizations 
has been well studied [3, EH EH , incidental benefits such 
as increased foraging efficiency have received little re- 
search attention. In contrast, research on foraging ef- 
ficiency has focused largely on independent individuals 
13h18|. and comparing foraging behavior across species 
19(. Therefore, to date, very few models have examined 
the potential effect that long-distance calling can have on 
movement behaviors and population distributions, and 
many open questions remain, particularly in the inter- 
relation between communication and optimal search for 
resources. For example, can this communication really be 
considered 'dependent' on the distribution of resources? 
How can communication facilitate group formation and 
identification of forage resources? Does a communication 
range exist that optimizes foraging efficiency? To what 
degree does search efficiency depend on the communica- 



tion mechanism? Finally, how does communication affect 
individual space use in a heterogeneous environment? 

In this letter, we address these questions with a gen- 
eral model of random search with two main ingredients: 
resource gradients and long-range communication. We 
first introduce the theoretical model and show how search 
time changes when foragers share information. Then 
we apply a specialized version of the model to the par- 
ticular case of acoustic communication among Mongo- 
lian gazelle, the dominant wild herbivore in the Eastern 
steppes of Mongolia, which must find each other and rel- 
atively small areas of good forage in a vast landscape, 
where noise travels substantial distances ■ We aim to 
explore whether acoustic communication in the Steppe 
could lead to the formation of observed large aggrega- 
tions of animals 2l|, and how search efficiency depends 
on the distances over which calls can be perceived. We 
wonder if the call length-scales that optimize search in 
real landscapes are biologically and physically plausible. 
To do this, we couple an individual-based representation 
of our model with remotely-sensed data on resource qual- 
ity in the Eastern Steppe. 

We consider TV particles which undergo a 2 — D Brown- 
ian random walk biased by the gradients of the landscape 
quality (local information), and by the interaction among 
individuals through a communication mechanism that is 
activated when good resources are found, thus provid- 
ing information on habitat quality in other areas (non- 
local information) . The dynamics of any of the particles 
a = 1, N is 

r a (t) = B g Vg(r a ) + B c VS(r Q ) + Va (t), (1) 

where rj a (t) is a Gaussian white noise term characterized 
by (Va(t)) = 0, and <»fe(*hs(* / )) = 2DS aP S(t - t'), with 
D the diffusion coefficient. The term B g Vg(r a ) refers to 
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the local search, where g(r) is the environmental quality 
function (amount of grass, prey, etc..) and B g is the 
local search bias parameter. Bq^7 S(r a ) is the nonlocal 
search term, with Bq the nonlocal search bias parameter 
and S(r a ) is the available information function of the 
individual a. It represents the information arriving at 
the spatial position of the animal a as a result of the 
communication with the rest of the population. This 
term makes the individuals move along the gradients of 
the information received. This is a function of the sum of 
pairwise interactions between the individual a and each 
one of its conspecifics, 

5 Q (r) = F[ J2 A[g(r f3 )}V(r a ,^)) . (2) 

F is an arbitrary perception function that must be set in 
each application of the model, V(r a , rg) is the interaction 
between the receptor particle a depending on its position 
r a and the emiting particle fixed at r^, and ^4[<?(rg)] is 
the activation function (typically, a Heaviside function) 
that indicates that the individual at rg calls the others 
if it is in a good habitat. 

From the Langevin equation (til), and following the 
standard arguments presented in [22|, [23| it is possible to 
write an equation for the evolution of the number den- 
sity of individuals, p(r,t). This approach will allow us 
to fix the parameters of the problem having a better un- 
derstanding of the role they are playing in the dynamics 
through a dimensional analysis. However, in the case of 
ungulates we are going to study later, it is not very suit- 
able to describe a population as a continuum since the 
number of individuals is not very high and the typical dis- 
tances among them is also large. Neglecting fluctuations 
the continuum equation for the density is 

= BV 2 p(r, t) + B g \7 [p(r, t)Vg(r)} + 
+B C V [p(r, t)VF (J dv'p(r, t)A[g(r']V(r, r'))] . (3) 

The total population is N — J Q p(r, t)dr, and remains 
constant since Eq. ([3]) is conservative in agreement with 
the microscopic model, where we do not consider death 
and birth processes. This equation is quite similar to the 
one derived in [24( to study the transport of interacting 
particles on a substrate. 

The model can, in principle, be applied to many situa- 
tions where animals obtain information from conspecifics. 
As previously stated, we wish to explore how foraging 
times are affected when individuals share information, 
but our model could also be generalized to the case of 
predators which use prey's signals to locate them. For the 
general case, we consider an identity perception function, 
and a Gaussian-like interaction. Manipulating its typical 
range via the standard deviation, a, we ask how the typi- 
cal communication distance affects the average efficiency 
of individuals searching for some target in space (areas 



of high-quality forage). We give an answer in terms of 
spatial distributions of individuals at long times start- 
ing from a random initial condition, and the mean first 
arrival time, r, to the targets. 

We begin with Monte Carlo simulations of the 
individual-based dynamics in ([T]) using a theoretical land- 
scape quality function, g(r), consisting of three non- 
normalized Gaussian functions, to ensure that g(r) G 
[0, 1], centered at different spatial points, and using pe- 
riodic boundary conditions. The available information 
function of the individual a depending on its position 
will be 

* expf- (r V r / )2 N ) 
S(r a )= ^Ml ~2^2 "> ( 4 ) 

where, as mentioned before, A[g(r)] is a theta Heaviside 
function that activates the interaction when the quality 
is over a certain threshold K, A[g(r)] = Q(g(r) — k). 

We observe that the first arrival time (Fig. Q] (right)) 
may be optimized with a communication range parame- 
ter, er, of intermediate scale. The number of individuals 
from which a given animal receives a signal will typically 
increase with the interaction scale. When it is too small, 
the receiving individual obtains too little information and 
thus exhibits low search efficiency (Fig. [IJ. Similarly, 
interaction scales that are too large lead to the receiv- 
ing individual being overwhelmed with information from 
all directions, also resulting in inefficient search (Fig. [I]). 
Only intermediate communication scales supply the re- 
ceiving individual with an optimal amount information 
with which to efficiently locate the callers and the high- 
quality habitat areas they occupy. 

The same behavior is also shown by the macroscopic 
equation ([3]) (left panel in Figure [IJ . Now r is defined 
as the time that passes until half of the population has 
found a target, that is f,, >K p(r,t)dr > N/2. We have 
integrated the equation (|3]) in I — D system of length 
L. We have used a single Gaussian patch of resources 
centered at L/2 and periodic boundary conditions for a 
random initial condition, which is equivalent to the case 
of an infinite system with equidistant high quality areas. 
We have taken the calling bias as being much stronger 
than the resource bias to make the nonlocal mechanism 
much more important in the search process and thus eas- 
ier to see how the communication range parameter affects 
the search time. We should remark on the differences 
between the 2 — D individual-based and the 1 — D deter- 
ministic density equation description. These differences, 
coupled with the parameter choice (stronger bias in the 
density equation) and the size of the systems used in each 
case, explain the different observed time-scales in the left 
and right panels of Figure [TJ The distributions of indi- 
viduals in the long time limit, show that all the animals 
are in good habitats, i.e., in areas where the maxima of 
the g function occur (not shown). 
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FIG. 1: Search time using the macroscopic equation (left) 
with B g = 0.50, B c = 0.75, D = 0.05 and the individual 
based description (right) with B g = 0.50, B c = 50, D = 0.75. 
K — O.S5g m ax in both panels. 



Next, we present the application of the model to the 
Mongolian gazelle (Procapra gutturosa) using remote- 
sensing data quantifying habitat quality in the Eastern 
Steppe of Mongolia. The Normalized Difference Vegeta- 
tion Index (NDVI) is one of the most widely used veg- 
etation quality estimators that can be calculated from 
satellite imagery. It is a useful index widely used to pre- 
dict habitat quality in grazing ecosystems pjj . and it 
has been also applied to gazelle habitat associations in 
the Mongolian Steppe .26}. NDVI is characterized by 
the function ffd(r), a continuous function taking values 
between (no vegetation) and 1 (fully vegetated). As 
the vegetation at low NDVI is too sparse, and at high 
NDVI it is too mature and indigestible, gazelles typically 
seek forage patches that fall intermediate NDVI values 
[26| . To make gradients of resources drive the movement 
of the individuals to regions with intermediate NDVI val- 
ues, we apply to the data the following linear transforma- 
tion: g(r) = g d {r) if g d (r) < 0.5, and g(r) = 1 - g d (r) if 
<7d(r) > 0.5. It defines a resources landscape with values 
between (0, 0.5) where represents both fully vegetated 
and no vegetation (i.e., low quality forage). 

We study an area of ~ 23000 km 2 in Mongolia's East- 
ern steppe [32] and assume that the resources remain 
constant in time during foraging. As the data were sam- 
pled with 500 m resolution, we have applied the resource 
function, obtained from the linear transformation of the 
NDVI data, g(r), to the nodes of a 2 — D regular lattice 
with 500 m spacing and reflecting boundary conditions. 
It is crucial now to properly choose the perception func- 
tion in order to realistically model the case of gazelles per- 
forming acoustic communication. It is well known that 
the sensitivity of the response of the ear does not follow a 
linear scale, but approximately a logarithmic one. That 
is why the bel and the decibel are quite suitable to de- 
scribe the acoustic perception of a listener. Therefore we 
choose an acoustic perception function of the form 



S(r a ) = 10 log 
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a A[g(rp)]V(r a ,rpy 
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where the sound calling of j3, V(r,rp), plays the role 
of a two body interaction potential, Iq is the low per- 
ception threshold (we take the value of a human ear, 



mals [27jj , and in any case, is just a reference value on 
which our results will not depend). The interaction po- 
tential mimicking acoustic communication is 
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considering that sound from an acoustic source attenu- 
ates in space mainly due to the atmospheric absorption 
(exponential term), and the spherical spreading of the 
intensity (47rr~ 2 contribution), and neglecting secondary 
effects [23| . Pq may be understood as the power of the 
sound at a distance of 1 m from the source. The ab- 
sorption coefficient, 7 is given by (Stoke's law of sound 
attenuation [27j |) 
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where r\ is the viscosity of the air, p the density, v the 
propagation velocity of the acoustic signal (which de- 
pends on the temperature and the humidity), and v 
is the sound frequency. We work under environmen- 
tal conditions of T = 20°C, and relative humidity of 
HR = 50%, which are quite close to the correspond- 
ing empirical values for the summer months from the 
Baruun-Urt (Mongolia) weather station, averaged over 
the last 4 years. These values give an absorption coeffi- 
cient of 7 « 10 °i/ 2 m _1 . The inverse of the absorption 
coefficient, 7 -1 , gives the typical length scale for the com- 
munication at each frequency, and thus plays the same 
role as the standard deviation, a, did in the Gaussian in- 
teraction used in the general model. From its functional 
dependence, different values of the frequency will modify 
the value of the absorption coefficient, and consequently, 
will lead to different communication ranges. Therefore, 
we will use sound frequency, v, as the control parameter 
of the interaction range. 

Using GPS data tracking the positions of 36 gazelles 
from June 2007 to August 2008, and assuming that they 
undergo Brownian motion in space, we estimate a dif- 
fusion constant of D — 74 km? day^ 1 . The bias pa- 
rameters in Eq. ([3]) have dimensions of [B g ] = m? s _1 
and [Be] — m 2 s , and to give empirically-based values 
to these parameters, we fix the total drift velocity as a 
constraint. Based on previous field work [26} we set 



Vdrift = B g Vg(v) + B c VS(r) = 10 km day' 



(8) 



The local search mechanism is responsible for short-range 
slow movements, while nonlocal communication gives 
rise to long and faster movements, and thus we require 
B g Vg(r) « 5 c VS(r). 

We have run simulations of a stochastic individual- 
based model with the dynamics of Eq. (JTJ for a one- 
month period for a population of N = 500 individu- 
als, which are initially distributed randomly on a 2 — D 
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FIG. 2: (Color online). Mean arrival time for 500 gazelles 
(averaged over 50 realizations). Parameter values: D = 



74 km 2 da\ 
13 km 2 day' 



B g = 2.6 x 10 _J fcm J day' 

= 0.70g m ax 



B c = 



regular lattice. We quantify the efficiency of the search 
for areas of high quality resources in terms of the mean 
first arrival time (averaged over 50 different initial condi- 
tions) of the population, and explore the dependence of 
this metric on the communication length, 7 -1 , or equiva- 
lently the frequency, v (Figure [5]). We obtain an optimal 
foraging time of 41 hours for 7 _1 of the order of kilo- 
meters (around 6 km). This result cannot currently be 
checked with data. However, switching to frequencies, 
the optimal search is obtained when gazelles communi- 
cate at a frequency of 1.25 KHz, which lies inside the 
measured interval of frequencies of the sounds emitted by 
gazelles, [0.4, 2.4] KHz [2o|,[2gj]. As in the general model, 
this means that the search is optimal when the receiving 
individual has an intermediate amount of information. 
Again, a lack of information leads to a slow, inefficient 
search, while an overabundance of information also de- 
creases foraging efficiency because the receiving animal 
is called from everywhere and gets lost in the landscape. 
These different regimes are also observed in the long time 
spatial distributions (i.e. efficiency of the search in terms 
of quality) shown in Figure [3] where we show the posi- 
tions of individuals after one month for different commu- 
nication lengths. At intermediate communication scales, 
v = 1 kHz, (Fig. [3] bottom left) all of the animals end 
up in regions with the best resources, regardless of where 
they started from. For smaller (Fig. [3] top) or larger 
(Fig. [3] bottom right) frequencies, some animals are still 
in low-quality areas at the end of the simulation period. 

In summary, our study clarifies some questions on the 
relationship between communication and optimal search 
for resources. The key result from our work is that, in 
general, intermediate communication distances optimize 
search efficiency in terms of time and quality. Individ- 
uals are able to find the best quality resource patches 
regardless of where they start from, opening new ques- 
tions about the distribution of individuals in heteroge- 
neous landscapes. The existence of maximum search effi- 
ciency at intermediate communication ranges is robust to 
the choice of functional form of V(r), allowing the model 
to be generalized to many different ways of sharing in- 
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FIG. 3: (Color online). Spatial distribution of 500 gazelles 
after 1 month, v = 0.1 kHz (top), v = 1 kHz (bottom left), 
v — 15.8 kHz (bottom right), where many individuals are 
in the same position. The size of the star is related to the 
number of animals. NDVI based resources landscape. 



formation. In the first part of the letter, we performed 
a theoretical study using a Gaussian interaction kernel 
and an identity perception function. In the second part, 
we focused on the case of acoustic communication among 
Mongolian gazelles, so that spreading of sound is mod- 
eled with a power law with an exponential cut off, and 
sound perception is a logarithmic function. We obtain 
an optimal calling distance empirically plausible, on the 
order of kilometers, similarly to the one observed in lions 
[30l | or hyenas 31 1 

intermediate 



Communication over intermediate scales results in 
faster search, and all the individuals form groups in areas 
of good resources. While this has obvious advantages in 
terms of group defense and predator swamping, it will 
also lead to rapid degradation of the forage (and thus 
habitat quality) at those locations. This is the problem of 
foraging influencing the patterns of vegetation, which will 
be treated in the future. Shorter-scale communication 
implies an almost individual search, which helps preserve 
local forage quality, but has clear disadvantages in terms 
of group defense strategies. On the other hand, longer 
scales lead to the formation of big groups (faster degra- 
dation of foraging), and animals need more time to join 
a group, which has negative consequences against preda- 
tion. Furthermore, acoustic communication scales signifi- 
cantly larger than the optimal scale for foraging efficiency 
identified here would be biologically implausible, even if 
ultimate group size (and not rate of group formation) was 
the most important aspect of an antipredation strategy. 
Exploring tradeoffs between group defense and individ- 
ual foraging efficiency in highly dynamic landscapes may 
open new research opportunities. 
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